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In this work, we model the zero-bias conductance for the four different DNA strands that were used in con- 
ductance measurement experiment [A. K. Mahapatro, K. J. Jeong, G. U. Lee and D. B. Janes, Nanotechnology, 
18, 195202 (2007)]. Our approach consists of three elements: (i) abinitio calculations of DNA, (ii) Green's 
function approach for transport calculations and (iii) the use of two parameters to determine the decoherence 
rates. We first study the role of the backbone. We find that the backbone can alter the coherent transmission sig- 
nificantly at some energy points by interacting with the bases, though the overall shape of the transmission stays 
similar for the two cases. More importantly, we find that the coherent electrical conductance is tremendously 
smaller than what the experiments measure. We consider DNA strands under a variety of different experimental 
conditions and show that even in the most ideal cases, the calculated coherent conductance is much smaller 
than the experimental conductance. To understand the reasons for this, we carefully look at the effect of deco- 
herence. By including decoherence, we show that our model can rationalize the measured conductance of the 
four strands, both qualitatively and quantitatively. We find that the effect of decoherence on G : C base pairs 
is crucial in getting agreement with the experiments. However, the decoherence on G : C base pairs alone 
does not explain the experimental conductance in strands containing a number of A : T base pairs. Including 
decoherence on A : T base pairs is also essential. By fitting the experimental trends and magnitudes in the 
conductance of the four different DNA molecules, we estimate for the first time that the deocherence rate is 6 
meV for G : C and 1.5 meV for A : T base pairs. 

PACS numbers: 87.14.gk, 73.63.-b, 87.15.A- 



I. INTRODUCTION 

Molecular electronic devices have been increasingly at- 
tracting attention from researchers working on both top 
down and bottom up approaches.^ ^ Among various molec- 
ular candidates, DNA has the attractive features of recog- 
nition and self-assembly."^ ^ Many efforts^ have been de- 
voted to understanding the mechanism of charge transfer 
and transport in DNA. From an electrical device perspec- 
tive, double barrier reso nant tunneling structures, spin spe- 
cific electron conductor,'^^^ field effect transistor,^^ trinary 
logic, "^optomechanical molecular motor,^^ negative differen- 
tial resistance,!^ detection of lesions by repair proteins^^ and 
dopinJE are in principle possible with DNA. Furthermore, 
recent work has also shown that it is possible to detect dis- 
eases by measuring the electrical conductivity of 
Researchers have also fou nd th at DNA can either be an 
insulator,^^ a semiconductor,^^^^ an ohmic conductor,^ ^ or 
a superconductor,^^ depending on a variety of conditions. 

It is broadly agreed upon that while the understanding of 
the electrical conductance of DNA has matured over the last 
two decades, explaining experiments or making predictions 
remains a challenge, compared to nanoengineered materials 
such as nanotubes and nanowires that have comparable di- 
mensions. Modeling of the conductivity of a DNA molecule 
in solution is complex because many factors contribute to the 
charge transport process. In wet DNA, water molecules are 
critical in influencing the molecular structure and hence the 
electronic properties. A-DNA usually has five to ten water 
molecules per base, while for B-DNA more than thirteen wa- 
ter molecules per base are preferred. At different hydration 
levels, different local densities of states (DOS) are obtained. 
Besides, because the backbone of a DNA molecule is neg- 



atively charged in its phosphate groups, it is believed that 
the molecule is surrounded by cations such as A/'a+, 
M^^+ and HsO^. Previous studies have shown that varia- 
tions in cation position^^ and type^^ can modulate the highest 
occupied molecular orbital (HOMO) and lowest unoccupied 
molecular orbital (LUMO) levels. Configuration changes (lat- 
tice vibrations/deformation) can also modify the conductivity 
of DNA molecules. Theoretical calculations have confirmed 
that changes in the distance and the angle between consec- 
utive bases can cause variability in conduction channels by 
influencing the hybridization between the tt orbitals of adja- 
cent bases. Additionally, the length and sequence play an 
important role in determining the conductivity of DNA. The 
sequence is important^^ because of the different ionization 
potential of the four bases, adenine (A), thymine (T), guanine 
(G) and cytosine (C) and sequence dependent conduction .E^SHl 

The conductivity of DNA molecules lying between two 
metal contacts has been studied by a few different groups. 
References 1 1713614 l[^50l have examined the conductance of 
DNA molecules whose coordinates are frozen based on the 
Landauer-Biittiker formalism. Their approaches can be clas- 
sified into two categ ories, the fi rst of which uses density func- 
tional theory (DFT).EISEIIil 

Because DFT calculation for 
long DNA strands is very time-consuming, researchers alter- 
natively first carry DFT calculations on short strands (typ- 
ically two to five base pairs), and extract the Hamiltonian 
and overlap matrices for the bases and their interactions with 
neighbors. One can then construct a larger matrix corre- 
sponding to the Hamiltonian for a longer strand from the sub- 
Hamiltonians obtained from the mentioned calculations on 
shorter strands. Another category is to use more simplified 
Hamiltonians that account for only one or two energy levels 
on each base and their interactions with energy levels in neigh- 
boring bases. 
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The modeUng of conductance in wet-DNA lying between 
metal contacts is a more difficult problem because the sur- 
rounding environment (water molecules, ions and conforma- 
tion) fluctuates with time. In this case recent work by Ref. 
[5l uses a combination of QM/MM to evaluate the role of 
fluctuating environment on DNA transport. These calcula- 
tions freeze the location of the atoms at each sampling point 
in the MD simulation and calculate the phase-coherent con- 
ductance using the Landauer-Buttiker approach. In addition 
to the above mentioned work that uses the Landauer-Buttiker 
approach, there has also been work that model DNA using a 
rate equation involving hopping and tunneling depending on 
the sequence involved.^ 

In this work, we study the conductance of dry DNA by care- 
fully investigating the role of decoherence. Decoherence in 
dry DNA molecule mainly arises from time-dependent fluc- 
tuations where the electrons will lose phase information to 
the environment. Electrons in dry DNA can lose phase in- 
formation by interacting with lattice vibrations and ambient 
electromagnetic fields. In this paper, we model t he de coher- 
ence using the phenomenological Biittiker probes. ^^^I^ We do 
this in the context of four different dry DNA strands consid- 



II. MODEL AND METHOD 



ered in the experiments of Ref. 40 and shown in Fig. 1(a) 



The four strands, which will be further discussed in section 
II, contain zero, one, three and five A : T base pairs in the 
middle part of a strand that otherwise consists of only G : C 
base pairs. We note that the method of including decoherence 
is not unique. Four perceptive recent papers include the ef- 
fect of decoherence in wet-DNA using completely different 
approaches.^^ "^^ ' ^^ ' ^^ ' Reference 57 represents the influence of 
atomic charges from the DNA backbone, water molecules and 
counterions with an empirical force field using a hybrid quan- 
tum mechanics-molecular mechanics (QM/MM) framework. 
Reference 58 uses a new statistical decoherence model devel- 
oped by their group. References l44l and l45l use the harmonic 
phonon bath to describe the dissipative environment. Refer- 
ence 46 includes the effect of energetic vibronic coupling via 
a full-fledged nonequilibrium Green's function approach. In 
our approach the DNA strands are fixed and charge transport 
in DNA is a decoherent process. A self-energy is used to rep- 
resent the decoherence. A similar approach has been adopted 
by Ref. ^59, with a simplified model Hamiltonian. In our cal- 
culation, we use DFT to get the full Hamiltonian which can 
describe the system more accurately. 

The remainder of this paper is organized as follows: In Sec. 
II we describe the method used for studying the conductance 
in detail. In this part, the Green's function formalism and 
Biittiker probes are reviewed. Then, we discuss the results 
in Sec. III. We first study the phase-coherent conductance of 
the four strands both with and without the backbone and find 
that the phase-coherent conductance is many orders of magni- 
tude smaller than the experimental values of the conductance. 
We then include decoherence in these strands to understand 
how decoherence changes the conductance. We show that by 
using two different decoherence rates one can provide reason- 
able agreement with experiments. Finally, we make a brief 
conclusion in Sec. IV. 



We consider four different DNA strands, in which the 
five base pairs in the center part are changed from G : 
C to A : T [Fig. l(a)| . The choice of our model 
system is motivated by the experiment in Ref. [40j 
which provides data for the room temperature conduc- 
tance of these four different strands under similar condi- 
tions, where the strands are connected to electrical con- 
tacts via thiol groups. The four strands we study consist of 
the following 15-base pair double- stranded DNA molecules: 
(a) GGCGCGCGGGCGGGC', (b) GGCGCGGAGGCG 
GGC- (c) GGCGCGAAAGCGGGC; (d) GGCGCAAA 
AACGGGC [Fig. |l(a)| . Conductances of these strands were 
measured in Ref. |40I under dry condition. The polycation 
spermidine was used to stabilize the double- stranded DNA 
molecules. In these systems, both the base pairing and base 
stacking of B-DNA is expected to remain unaffected. 

Our method involves the following five steps. First, we ob- 
tain the atomic coordinates for the double- stranded B-DNA 
using the Nucleic Acid Builder (NAB) software package.!^ 
To study the effect of the backbone, we consider two types 
of strands - the native ones with the backbones and the ones 
whose backbones are replaced with the hydrogen atoms. The 
sketch of the latter case is shown in Fig. |l(b)| where the 
arrows indicate the hydrogen atoms that replace the back- 
bone. Second, we use Gaussian to obtain the Hamil- 
tonian and overlap matrices, Hq and ^o, where the B3LYP 
functional and the 6-3 IG basis set^^^^ are adopted. We then 
transform the system to an orthogonal basis using Loewdin 
transformation,!^^^ 



H = U^HiU 



(1) 



(2) 



InEq. ([2|,[/ is a block diagonal matrix. To obtain U, we first 
diagonalize every diagonal sub-block of Hi and then arrange 
the eigenvectors in the order of DNA bases. While the diag- 
onal blocks of Hq are full matrices, the diagonal blocks of H 
are diagonal matrices and the dimension of the diagonal block 
matrices is equal to the number of orbitals used to represent 
that base. Physically, the diagonal blocks of H correspond to 
the localized energy levels of each DNA base and off-diagonal 
blocks correspond to interactions between different bases. If 
we look at the values of the off-diagonal elements, we find 
that the interaction corresponding to energy levels at the two 
nearby bases is large while that corresponding to two far away 
bases is relatively small. 

The transmission through the molecule is computed us- 
ing the Green's function approach,^^ with the Biittiker probes 
to account for the decoherence. We calculate the retarded 
Green's function, which is defined by. 



(3) 



where H is the Hamiltonian shown in Eq. ([2|. The self- 
energy of the left (right) contact ^l{r) represents the cou- 
pling of the DNA to left (right) contact through which charge 
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enters and leaves the DNA. The self-energy due to the phase- 
breaking Biittiker probes is . 

The self-energy due to the contacts is the third step of our 
calculations. The correct expressions for the self-energy due 
to the contacts are = TdlQlTld and = ToRgnTRD. 
Here ql and qr are the surface Green's functions of the 
contacts, Tld and Tjid are the coupling between the left 
and right contacts and the device, with Tdl = and 
Tdr = Tj^D' Computing these self-energies for realistic con- 
tacts is challenging. To obtain an accurate answer, the inclu- 
sion of contact surface atoms with optimized structure and 
their influence on the device molecule are required though 
some progress has been made recently in this area.^^ In this 
paper, we neglect the real part of the self-energy and further 
set ^L{R) = —i^L{R)/'^^ where ^l{r) is treated as an energy 
independent parameter. Mathematically, the coupling matrix 
is diagonal with the non-zero elements representing the cou- 
pling strength. This approxin iation h as been adopted by other 
researchers in DNA transport,'^^'^^'^ especially at small biases 
as long as the value of ^l{r) is close to experiments. We 
verify this approximation by trying different ^l{r) values. 

The fourth step involves the inclusion of Biittiker probes. 
Biittiker probes are fictitious probes which extract electrons 
from the device and re-inject them after phase breaking, as il- 
lustrated in Fig. [2] The net current at each Biittiker probe is 
zero. Similar to the effect of the left and right contacts, the 
effect of Biittiker probes is also included as a self-energy, E^, 
as shown in Eq. ([3|. Specifically, we use to represent the 
decoherence at the ith probe, = XI ^i- In our model, 

i 

is controlled by an energy-independent coupling strength be- 
tween the probe and the coherent system F^, = —iVi/2. In 
the calculations, we attach the Biittiker probes to the energy 
levels of the bases that correspond to the diagonal elements of 
H. The strength of coupling of electrons from the DNA to the 
Biittiker probe is assumed to depend only on whether it is a 
G \ C ov A : T base pair. Note that at first we include the 
Biittiker probes only in the^ : C base pairs as A : T base 
pairs are tunneling barriers.!^ We find that only including de- 
coherence on G : C is not enough to explain the experiment. 
The decoherence on A : T also plays a role in conduction, 
indicating that A : T is not a coherent static barrier. 

The fifth step uses the Green's function approach to cal- 
culate the conductance for various coupling strengths of the 
Biittiker probes. The structure of the DNA is assumed to be 
static in the conventional B form throughout the calculation. 
The effect of fluctuations in the lattice and environment is in- 
cluded only via the Buttiker probes. In the following part of 
this section, we summarize the steps of obtaining the effec- 
tive transmission and conductance after including the effect 
of Biittiker probes using the D' Amato-Pastawski model.^^ 

Assuming that the total number of DNA bases is N (here, 
N = 30), the number of Biittiker probes is A^5 = TV - 2. At 
the low-bias region, the current at the i th probe is, 



li = -^^Tijl^i - Ni^ = 1,2, • • . ,iV (4) 
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FIG. 1: (Color online) (a) Sketch of the four DNA sequences: Seq. 
1, Seq. 2, Seq. 3 and Seq. 4, where each base is numbered. Yel- 
low regions indicate the ends connecting to the contacts, (b) Atomic 
structures of the hydrogen atoms terminated DNA bases. Arrows 
indicate the hydrogen atoms that replace the backbones. 



The factor 2 is used to account for the spin degeneracy. In 
the above equation, Tij is the transmission between the zth 
and jth probes calculated by T^j = F^G^F^G^, where 
is the advanced Green's function, G^ = (G^)^ Because the 
transmission coefficients are reciprocal, we have Tij = TjiW^ 




FIG. 2: (Color online) Sketch of the device with Biittiker probes. 
Yellow regions indicate the ends connecting to the contacts. Blue 
regions indicate DNA bases and black dots are Buttiker probes. Red 
arrows describe the behavior of electron at Biittiker probes: first ex- 
tracted from the device and then re-injected back. 

Using the condition that the net current is zero at each 
Biittiker probe, Eq. Q gives us Nij independent formulas, 
from which we can express the chemical potential of the ith 
Biittiker probe, /i^, 

Nb 

Mi-ML = E^7'T,i,](/i^-/i^),i = l,2,--- (5) 

In the above equation, is the inverse matrix of W, 

whose elements are given by,!^ W^j = [(1 — Rii)Sij — Tij{l — 
Sij)], where Ra is the reflection probability at probe i, which 
is given by Ra = 1 — Y^^^i Tij . The currents at the left 
contact II and right contact Ir are not zero. Because of the 
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conservation of the electron number, II -\- Ir = 0. Express- 
ing the current at the left contact in terms of the difference 
between the chemical potentials at the two contacts fiL and 



(6) 



Comparing Eq. (|6| and the expression for the current at 
the left contact given by Eq. Q, we can write the effective 
transmission between left contact and right contact as. 



(7) 



The first term Tlr is the coherent transmission from left 
contact to right contact while the second term describes the 
effect of decoherence. From Eq. ([6|, the zero-bias conduc- 
tance is approximately. 



(8) 



Go = ^ ^ 7.75 X 10 ^ is the quantum conductance. 



III. RESULTS AND DISCUSSION 

The sketch of the four 15 -base pair DNA sequences that we 
model is shown in Fig. 1(a) In the experiment, thiol groups 
are used to connect the 3' end of the DNA molecules to the 
metal contacts, because thiol groups are expected to provide 
a strong coupling between the DNA molecule and the metal 



contact. In this work, we use two broadening matrices, Tl 
and Tji to represent that coupling. We vary F l and F^^ be- 
tween 50 and 100 meV, and find that there is no significant 
change in our results. For the results presented here, we set 
F^ = Fi^ =100meV. 

In this section, we analyze the effect of two factors in 
DNA charge transport using the Landauer-Biittiker formal- 
ism, namely, the backbone and the decoherence. In Section 
III A, we compare the transmission for the four strands with 
and without the backbones in the phase-coherent limit. We 
find that the backbone does not change the transmission qual- 
itatively. However, it can affect the transport significantly at 
some energy points by interacting with the atoms in bases. In 
Section III B, we discuss the role of decoherence. We observe 
that the coherent transmission is too small to explain the ex- 
periment. We infer that the low value of conductance can be 
due to a number of reasons: poor coupling to contacts, po- 
sition of Fermi energy or incorrect bandgap. We will show 
below that even for the best possible values of these parame- 
ters, the phase-coherent conductance is smaller than the exper- 
imental conductance. We will then show that by including ap- 
propriately chosen decoherence rates the calculated conduc- 
tance values are comparable to the experiment. The strength 



of decoherence rate is difficult to compute abinitio. Here we 
have considered a variety of values to lend some insight into 
the ones which may be experimentally feasible. In the process 
of modeling decoherence, we have also found that it is possi- 
ble to explain the experiments only if we assume that there is 
also decoherence in the A : T barriers. That is, the A : T 
barriers cannot be assumed to be static rectangular barriers. 



A. The Role of the Backbone 

To study the effect of backbone, we have calculated trans- 
mission for two cases - DNA strands with and without back- 
bones in the phase-coherent limit. For the strands with back- 
bone, we use the positively charged sodium ions Na~^ to neu- 
tralize the backbone. To determine the geometry of the DNA 
strand after being neutralized, we first carry DFT calculation 
on a short strand containing five G stacking bases with sodium 
ions placed nearby the phosphate group by setting the initial 
distances based on Ref. |42| The atom coordinates for DNA 
molecule are fixed while the positions of Na~^ ions are re- 
laxed. We approximate the optimized position of the sodium 
ion around the backbone of the middle G in this short strand 
to be the positions of the thirty sodium ions in the long 15- 
base pair strand. We then carry out DFT calculation on the 
15-base pair strand to obtain the self-consistent Hamiltonian 
and overlap matrices. For the strands without backbone, we 
simply delete the backbone and then terminate the base with 
hydrogen atoms using GaussView 5^ 

The transmission of the four DNA strands is calculated in 
the phase-coherent limit as shown in Fig. [3] For the strands 
with backbone, we model the contacts via two approaches - 
charge is injected and extracted at the: (i) base and backbone 
(blue solid) and (ii) base only (magenta dash). Transmission 
for the strands without the backbone is shown by the red dash- 
dot curve of Fig. |3] We note that the overall shapes of trans- 
mission with and without the backbone qualitatively seem to 
be similar but the two curves are shifted in the energy axis. To 
obtain a quantitative understanding, we shift the energy axis 
of the transmission for the strand without the backbone to the 
right-hand side by 0.85 eV for Seq. 1 and 2, and 1 eV for Seq. 
3 and 4, as shown in Fig. [4] A close inspection of Fig. [4] re- 
veals that the difference in magnitude of transmission for the 
strands with and without the backbones can be up to 5 orders 
for Seq. 1, 6 orders for Seq. 2, 7 orders for Seq. 3 and 8 orders 
for Seq. 4 at some energy points. 

We study the effect of backbone further by plotting the 
HOMO orbitals for the four strands. HOMO orbitals for 
DNA strands with the backbones have been studied by a few 
groups. ^^^^^"^ It has been shown that the spatial distribution 
of HOMO orbitals depends on a variety of conditions, such as 
the position of counterions relative to the phosphate groups^ 
and the hydration levels. In addition, Ref. |^ states that 
the HOMO orbitals can temporarily have a large weight on 
the backbones as a function of time. In our calculation, we 
find that the HOMO orbital is distributed over both the bases 
and the backbone, as shown in Fig. [5] We conclude from 
the above observations related to transmission ( Fig. |4]) and 
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FIG. 3: (Color online) Transmission vs energy for the four DNA 
strands with and without the backbones in the phase-coherent limit. 
The transmission for the strands with the backbones is calculated by 
two methods- charge is injected and extracted at the: (i) base and 
backbone (blue solid) and (ii) base only (magenta dash). Transmis- 
sion for the strands without the backbone is shown by the red dash- 
dot curve. The transmission through the base only is slightly smaller 
than that through the base and backbone and the overall shapes of 
transmission with and without the backbone qualitatively seem to be 
similar but the two curves are shifted in the energy axis. 
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keeping the overall shape similar. We also note that the co- 
herent transmission through the HOMO-LUMO gap also de- 
pends on the inclusion of the backbone. This finding indicates 
that when modeling the coherent charge transport in DNA 
molecules, one should take the effect of backbone into ac- 
count. In the remainder of this paper, we will focus on the 
strands with the backbones, injecting and extracting from and 
into the strands via both backbone and base. 
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FIG. 5: (Color online) Isosurfaces of HOMO orbitals for Seq. 1 (a), 
Seq. 2 (b), Seq. 3 (c) and Seq. 4 (d). The red color represents the 
positive part of the wave function and the green color represents the 
negative part of the wave function. The HOMO orbital is distributed 
over both the bases and the backbone. 



FIG. 4: (Color online) Transmission vs energy for the four DNA 
strands with and without backbone in the phase-coherent limit. The 
transmission for the strands with backbone is shown in blue curves, 
while that for strands without backbone is shown in red with the en- 
ergy axis shifted 0.85 eV for Seq. 1 and 2, and 1 eV for Seq. 3 and 
4. The difference in magnitude of transmission for the strands with 
and without the backbones can be significant at some energy points. 



wave function (Fig. |5]) that while most of the contribution to 
charge transport comes from the bases, the phosphate groups, 
sugar rings and sodium ions on the backbone alter the elec- 
tronic structure sensitively by interacting with the bases. This 
leads to mismatches in the coherent transmission for strands 
with and without the backbone at some energy points while 



B. The Role of Decoherence 

Besides the interaction between the molecule and the metal 
contacts, the position of the Fermi level Ej relative to the 
molecular energy levels is also important in understanding the 
flow of charge. Ef depends sensitively on the surface con- 
dition and work function difference between two materials. 
Experimentally, it is very difficult to determine the position of 
Ej. Here, we explore the conductance when the Fermi level 
is in the HOMO vicinity. At the low-bias region, only elec- 
trons with energies nearby the Fermi level contribute to the 
transport. Thus, the transmission T{E) = T{Ef). From DFT 
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calculations, the HOMO levels of the 4 strands have been de- 
termined to be -3.18 eV, -2.84 eV, -2.91 eV, -2.85 eV, respec- 
tively. 

The zero-bias conductance versus the position of Fermi en- 
ergy is shown in Fig. |6] Comparing the conductances of Seq. 
2, Seq. 3 and Seq. 4, we find that the A : T base pair plays 
the role of a barrier. This is consistent with the fact that in our 
calculation the ionization potential (IP) of an isolated G : C 
base pair is 7.12 eV while that of the A : T base pair is 7.68 
eV. The similarity of the conductance for Seq. 1 and Seq. 2 
can be attributed to the fact that Seq. 1 has a G : C base 
pair in the sixth position and a C : G base pair in the seventh 
position which can cause interstrand hopping, and results in 
lower transmission while in the Seq. 2 both the sixth and sev- 
enth positions are placed with G : C base pairs. However, the 
presence of A\T base pair in the eighth position cancels out 
this advantage. More importantly, in this coherent model, we 
find that the conductivity is much smaller than the experimen- 
tal results. At the low-bias region, the experimental value^^ 
of the conductance for Seq. 1, Seq. 2, Seq. 3 and Seq. 4 is 
around 5 x IQ-^^ S, 2 x IQ-^^ S, 3.5 x 10"^^ S and 6 x IQ-^^ 
S, respectively. However, in the coherent case of our calcula- 
tion, the conductance is only about 10 "^^iS for Seq. 1 and 
10"^^ S for Seq. 2, Seq. 3 and Seq. 4 when the Fermi level 
is nearby the HOMO. For Seq. 1, the coherent result is about 
10^^ times smaller than the experiment; for Seq. 2, 10^^ times 
smaller; for Seq. 3, 10^^ times smaller and for Seq. 4, 10^^ 
times smaller. 
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FIG. 6: (Color online) Conductance vs Fermi energy for the four 
DNA strands without any decoherence. The coherent conductance 
is orders of magnitude smaller than experiment irrespective of the 
locations of Fermi levels for all the four strands. 

The huge difference between the coherent conductance and 
experimental values indicates that there are other mechanisms 
that play a role in determining the conductance. We rule out 
the location of Fermi energy because irrespective of where the 
Fermi energy is in Fig. [6] (including well within the HOMO 
band), the conductance is significantly smaller than in exper- 
iments. We also rule out the nature of the coupling to con- 
tacts because we have verified that our results change by small 
amounts when the coupling to the contacts F l and F^^ is var- 
ied from 50 to 100 meV. Another reason for the larger experi- 
mental conductance could be the possibility of having a large 



number of strands, which is extremely unlikely in the break 
junction geometry of Ref. 40 

The mismatch in energy level between neighboring bases 
combined with the small interbase coupling leads to the 
smaller than experimental conductances in the phase coherent 
calculations. Decoherence due to interaction with the envi- 
ronment broadens these energy levels and can lead to a larger 
conductance. To account for decoherence, we use Btittiker 
probes as presented previously. 

As shown in Fig. [6j for hole transport, G : C base pair 
is preferred while A : T base pair is considered as a barrier 
because the transmission decreases with increasing number of 
A : T base pairs. As the transit time through a rigid barrier 
A : T is small, we first include decoherence only on the G : C 
base pairs with the coupling strength F^ = 5 meV, 6 meV and 
10 meV, as shown in Fig. [7] Compared with the ballistic 
transport, one can find that after including Btittiker probes the 
conductance is smoother for the four strands. The small peaks 
in the ballistic conductance are smeared out due to the broad- 
ening of energy levels. For Seq. 1, we find that when we add 
5 meV decoherence only on the G : C base pairs, the con- 
ductance increases by around 10^^ times around their HOMO 
levels, while for Seq.2, Seq. 3 and Seq. 4, the increase is 10^^, 
10^^ and 10^^ respectively. The tremendous enhancement in 
the conductance after including Btittiker probes suggests the 
crucial role of decoherence in charge transport through DNA 
molecules. 

In addition, we find that for strands containing A : T base 
pairs, the same values of decoherence rates can enhance the 
conductance more effectively for Seq. 2 than Seq. 3 and Seq. 

4. We observe that when the decoherence on G : G base pairs 
is 5 meV, as the number of A -.T base pairs increases from 1 to 

5, the increase in conductance decreases from around 10^^ to 
10^^. This is because in Seq. 3 and Seq. 4, where the numbers 
of A : T are 3 and 5, respectively, the wider barriers begins 
to play a larger role than in Seq. 2. We also find that when 
the decoherence on G : G base pairs changes from 5 meV to 
10 meV, the conductance increases around 4 times for Seq. 1 
and Seq. 2 around the HOMO levels. However, for Seq. 3 
and Seq. 4, when the decoherence value on G : G varies from 
5 meV to 10 meV, the conductance only increases 1.5 times. 
More importantly, for Seq. 1 and Seq. 2, as the Btittiker probe 
coupling strength changes, we can fit the experimental con- 
ductance by setting the Fermi energy to be a particular value. 
For instance, when the decoherence on G : G is 5 meV if we 
set the Fermi level to be -3.06 eV for Seq. 1 and -2.82 eV for 
Seq. 2, the computed conductance is comparable to that of 
experiment. However, for Seq. 3 and Seq. 4, where the wide 
A : T barriers are present, irrespective of the strength of deco- 
herence chosen (varying from to 10 meV), the conductances 
are still too small when compared with the experiments. For 
example, for Seq. 3 the best conductance value is 10 times 
smaller than experiment and for Seq. 4, the best value is 100 
times smaller. This indicates that the decoherence on A : T 
barriers is also important, that is, the A : T barrier is not a 
static barrier for hole transport. For this reason, we study the 
effect of decoherence on A : T base pairs to get better agree- 
ment with the experiments. 
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Fig. [8] shows the conductance versus Fermi energy with 
different decoherence values on both G : C and A : T base 
pairs. It is found that for Seq. 2, the decoherence on G : C 
base pairs is more important than that on A : T base pair. 
For example, when the decoherence on A : T base pair is 
fixed as 2 meV, the conductance increases 1.5 times with the 
decoherence on G : C base pairs increasing from 5 meV to 
6 meV. However, if the decoherence on G \ C base pairs is 
fixed as 6 meV, the coductance stays unchanged at the HOMO 
level when the decoherence on ^4 : T base pair is varied from 
1 meV to 2 meV. In contrast, for Seq. 3 and Seq. 4, the 
decoherence on A : T base pairs are much more important. 
Take Seq. 3 for example. If the decoherence on A \ T base 
pairs is fixed as 2 meV and the decoherence on G : C is varied 
from 5 meV to 6 meV, the conductance only changes around 
1.1 times. Similarly, if the decoherence on A : T is fixed as 1 
meV and the decoherence on G : C is changed from 6 meV to 
10 meV, the conductance also only changes around 1.1 times. 
However, if the decoherence on G : G is fixed as 6 meV and 
the docoherence on A : T is changed from 1 meV to 2 meV, 
the conductance increases around 2.8 times. 

We find that for the four strands a decoherence strength of 
6 meV on G : G base pairs and 1.5 meV on A : T base pairs 
give a relatively good agreement with experiments. The re- 
sults are summarized in Fig. [9] where the experimental values 
of conductance from Ref . |40| are also shown, with horizon- 
tal lines. We find that for all the four strands, the computed 
conductance lies in the ball park of the experimental values. 
That is, the qualitative trend of the conductance variation for 
the four sequences matches well with the experiments (Con- 
ductance decreases from Seq. 1 to Seq. 4). More importantly, 
in contrast to the phase coherent results in Fig. |6j where the 
conductance values are orders of magnitude smaller than the 
experiments, the conductance values in Fig. |9]are now com- 
parable to the experiments. This qualitative match holds well 
irrespective of the Fermi energy as can be seen in Fig. |9] Al- 
though for Seq. 4, the computed conductance is about 2 times 
larger than the experimental result around the HOMO level, 
we still think that the fitting is good because of the uncertain 
factors in the experiments. 

These results indicate that the decoherence is important in 
experiments involving dry DNA, and this might be one of the 
reasons why the theory has difficulty in explaining experi- 
ments. 

Recently, there have been proposals to calculate the con- 
ductance of DNA by calculating the phase coherent conduc- 
tance for many points (each of which provides a static set of 
distinct coordinates) along the trajectory of and MD simula- 
tion and averaging the coherently calculated conductance.^ 
These coordinates include the effect of time varying ions and 
water molecules, apart from thermally induced lattice vibra- 
tions. It would be interesting to compare results from the 
two methods for the case of dry DNA to see if one can match 



experiments on the linear response conductance and current- 
voltage characteristics. To start with one could compare the 
average DOS obtained by the two methods. Also, work in 
Ref. 5 1 calculates the phase coherent conductance of DNA 
for many realizations of coordinates obtained from MD cal- 
culations of DNA in water. Such a process may give rise to an 
effective broadening of energy levels. 



IV. CONCLUSION 

In summary, we have modeled the zero-bias conductance 
of dry DNA using a combination of density functional the- 
ory and the phenomenological Biittiker probes to account for 
decoherence. We first explore the effect of the backbone in 
charge transport by comparing the coherent transmission for 
strands with backbones and strands whose backbones have 
been deleted. We find that the DNA backbone can affect 
the coherent transmission significantly at some energy points. 
However, the overall shapes of the transmission are similar for 
the strands with backbones and strands without backbones. 
More interestingly, we find that the calculated conductance 
using phase coherent transport is orders of magnitude smaller 
than experiments, and that including the effect of decoherence 
on G : G is crucial because it broadens the energy levels of 
DNA and significantly enhances the conductance. By compar- 
ing with experiments, we find that quite a large decoherence 
is required even on the A : T base pairs, even though they be- 
have as barriers in the energy ranges of importance. It is also 
worth noting that for the G : G rich strands (strands which 
contain one A : T base pair), the decoherence on G : G base 
pairs is more important than that on A : T base pairs, while 
for the strands contain three or five A : T base pairs, the deco- 
herence on A :T plays a more significant role. By analyzing 
the four different DNA strands, it has been determined that 
the decoherence strength due to Biittiker probes is approxi- 
mately 6 meV for G : G base pairs and 1.5 meV for A : T 
base pairs. While the phenomenological Biittiker probes is 
able to explain the experiments qualitatively at low biases, a 
more accurate description of vibronic coupling would be nec- 
essary to explain experiments at large biases. This is a signifi- 
cantly more difficult problem requiring knowledge of both the 
vibronic modes and the coupling strength in DNA, which is 
well beyond the extent of this work. 
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FIG. 7: (Color online) The conductance versus Fermi energy with 5 
meV, 6 meV and 10 meV decoherence only on G : C base pairs for 
(a) Seq. 1, (b) Seq. 2, (c) Seq. 3 and (d) Seq. 4. The conductance 
values for Seq. 3 and Seq. 4 around their HOMO levels are still 
too small to explain the experiment, suggesting the decoherence on 
A : T barrier is also important. 
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FIG. 8: (Color online) The conductance versus Fermi energy with 
different decoherence values on both G : C and A : T base pairs for 
(a) Seq. 2, (b) Seq. 3 and (c) Seq. 4. The decoherence on A \ T base 
pairs can increase the conductance effectively, especially for Seq. 3 
and Seq. 4. 
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FIG. 9: (Color online) The conductance versus Fermi energy with 
6 meV decoherence on C : C and 1.5 meV decoherence on A : T. 
The conductance has been increased effectively. The conductance 
values for the 4 strands around their HOMO levels are quantatively 
comparable to the experimental results, demonstrating the impor- 
tance of adding decoherence on both G : C base pairs and A : T 
barrier. 



